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We present preliminary results of an excited state spectroscopy calculation in the 2-d lattice Gross- 
Neveu model. We address the construction of suitable interpolators for the variational method and 
their overlap with excitations. We comment on the role of the eigenvectors as a tool for matching 
scattering states on lattices with different volumes. 
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1. Introduction 

Lattice QCD is by now a mature field which produces reliable quantitative results from ab-initio 
calculations. With refined techniques the range of quantities that can be computed on the lattice is 
increasing continuously. 

Observables that recently have started to see a lot of attention in the lattice community are 
masses of excited states of light hadrons (see e.g. Excitations, however, are quantities that are 
notoriously difficult on the lattice. The reason is that the energies W n of excited states appear only 
as sub-leading exponentials in the spectral decomposition of Euclidean 2-point functions, 

Qj(t) = (Oi(t)Oj(0f) = £(0|O,-|n)(n|O]|0) e -^" . (1.1) 

n 

Here Oj, Oj are operators with the quantum numbers of the state one is interested in. The sum runs 
over all physical states \n) with these quantum numbers, and W n are the corresponding energies 
(W\ < W2 < W3 ...). It is obvious, that the exponential decay of the 2-point function is dominated 
by the energy of the ground state W\ , while the energies of the excited states ^2,^3... appear in sub- 
leading exponentials and thus their contribution to the 2-point function is exponentially suppressed 
relative to the ground state. 

Although other proposals for the extraction of the excites state energies exist [Qj, the most 
commonly applied technique is the variational method []3j, Q]. Here one uses a large basis of inter- 
polators Oj, i = 1,2, ...N and considers the generalized eigenvalue problem of the N x N correlation 
matrix C(t), 

C(f)vW = A(0 W C(/ )vW . (1.2) 
The ordered eigenvalues, A' 1 ) > A^ 2 ' > A' 3 ) > behave as [Q] 

AW(f) = e -^ w " [\ + @{e- (t - to)K )] , (1.3) 

where W n is the energy of the n-th state and A„ the distance of W n to the neighboring energy level. 
Effectively the ground and excited states are disentangled, and each energy appears in an individual 
eigenvalue. The variational method can only be as good as the basis of interpolators Oj one uses 
and different construction principles were proposed [[[]]. 

Even with the variational method, the calculation of excited states is a challenging enterprise. 
Several non-trivial steps are involved: One has to construct a suitable basis of interpolators and 
select from those a subset which optimizes the quality of the signal, i.e., provides the cleanest 
effective mass plateaus. For an extrapolation to the infinite volume, continuum and chiral limits, the 
signals from several lattices with different size, coupling and mass have to be matched. Here it turns 
out that the eigenvectors v W provide an important marker for identifying the states. Finally, in a 
fully dynamical simulation hadrons may decay. This leads to coexistence of bound- and scattering 



states which have similar energies W„ and thus mix in the spectral sums (|1.1|). Because scattering 
states have relative momentum, they can be identified by a finite volume analysis, since the minimal 
possible momentum is inverse proportional to the spatial extent of the lattice 

Although the methods for solving the technical problems of an excited spectroscopy calcula- 
tion are known, so far the typical analysis produced only the masses for one or two excitations and 
most papers are still in the quenched approximation, where hadrons cannot decay. 
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There is a long standing tradition of studying problems of theoretical physics in low dimen- 
sional quantum field theories. This allows to test new concepts and techniques in an environment 
where often the situation is simpler and, when the lattice is used, large statistics is easy to generate 
(examples for such 2-d studies are [Q], [^]). The experience gained in 2-d is then often an important 
guideline for the more demanding case of full QCD. 

In this paper we follow this tradition and present first results of a lattice study of excited states 
in the 2-d Gross-Neveu model [^|. This model is particularly attractive for excited state spec- 
troscopy, since in the infinite flavor limit, analytic results suggest that the masses of the excitations 
are integer multiples of the ground state mass. Our project aims at extracting with the variational 
method as many excitations as possible and to study their volume dependence to learn about scat- 
tering. 



2. Lattice simulation of the Gross-Neveu model 
2.1 Setting of the calculation 

We use the Wilson action for discretizing the fermions on a 2-d lattice A. The total action is 

x,yeA xeA 

The fields \j7, y are vectors of Nf two-component spinors, one for each flavor, and we use ma- 
trix/vector notation for all indices. <I> denotes a one-component real scalar field, which generates 
the 4-fermi interaction of the model through a Hubbard-Stratonovich transformation. The Dirac 
matrix D is diagonal in flavor space and reads 

D(x,y) = [ m + y/g *(„) ] % - f iZ^ . (2.2) 

n=±i 1 

The matrices are given by the Pauli matrices, /i = di, 72 = 02, 7s = 03. 

The dynamical simulation of the model is done with standard Hybrid Monte Carlo methods 
and we work with Nf = 2,4 and 6 flavors. Our L\ x L2 lattices have sizes L\ = 10. ..24, L2 = 
32 ... 64 and we typically use four values of the coupling constant g and several different masses 
m. Statistical errors are determined with the Jackknife-method. We find that statistics of a few 
hundred configurations per ensemble are sufficient for statistical errors in the one percent range. 
Throughout this paper we set the lattice constant equal to one and all results are in lattice units. 
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2.2 Interpolators 

For our analysis of excited states we use the following set of flavor singlet interpolators: 
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(2.3) 



To project to vanishing total momentum, we sum over the spatial index x, a step which we suppress 
in the list of Eq. (23) for notational convenience. The interpolators are built from field variables 
at different lattice sites, and the integer valued parameters n and m determine how many steps in 
spatial direction the field variables are displaced. 

Furthermore it is possible to have a relative minus sign between the displaced fields which 
gives rise to a derivative type of fermion source. In order to have the same quantum numbers as the 
combination with a plus sign, an additional factor of 71 appears. In 2-d all products of 7-matrices 
can be simplified to a single 7^, or 1 due to the algebra of the Pauli matrices. 

Our interpolators have negative parity, i.e., they acquire a minus sign under the parity transfor- 
mation 

y(x,t) y(x, P = 72 V(-* 5 > Wfa P = V(~x,t) 72 • (2-4) 
Furthermore all interpolators are eigenstates with C = — 1 of the charge conjugation 

Y( x ,t) Y(x,t) c = C 1 \j/{x,t) T , \j7(x,t) -^\j7(x,t) c = -\j/(x,t) T C , (2.5) 

where C is the charge conjugation matrix obeying Cy^C^ 1 = — J^, which can be chosen as C = iji. 
We stress that in 2-d no angular momentum exists and thus P and C completely characterize our 
flavor singlet states. 

In our analysis we use two different choices for the displacement parameters n and m. We set 
n = 3 and m = 3 to obtain the 9 interpolators 0\ ... Ow as listed in ( |2.3| ) (note that there is no 
interpolator O7). A second choice, n = 4 and m = 2, gives rise to our interpolators On ... O20 
labeled accordingly. Thus we can use correlation matrices with a maximum size of 16 x 16. 

The implementation of the correlation matrix was checked in two independent programs and 
compared to the results for the free case computed from Fourier transformation. Within error bars, 
the correlation matrix was found to be hermitian as expected. The decay properties of the individual 
entries in the correlation matrix are either of the cosh- or sinh-type. 
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Figure 1: The diagonal entries C,, (f ) of the correlation matrix (l.h.s. plot) and the corresponding effective 
energies (r.h.s) for a 16 x 64 lattice with g = 0.05, m — 0.05 and Nf = 2. The numbers in the legend are the 
labels of our interpolators. 



3. Results 

For a first assessment of the data we consider plots of the correlators and of effective energies, a 
term which is more suitable than the usual "effective masses", since some of our states will turn out 
to be scattering states, such that their energy does not correspond to a single mass. The effective 
energies are defined as 

H-«(, + 1 /2)=,„^, (3.1) 

where c(t) is either an individual entry of the correlation matrix, or one of the eigenvalues of the 
generalized eigenvalue problem ( |1.2[ ). 

In the l.h.s. plot of Fig. [I] we show the diagonal entries C„ (f) of the correlation matrix on a 
logarithmic scale. The r.h.s. plot contains the corresponding effective energies. The data are from 
16 x 64 lattices at g = 0.05, m = 0.05 and N f = 2. 

It is obvious from the l.h.s. plot that on the logarithmic scale the correlators C„ (f) fall on a 
straight line beyond t ~ 3, i.e., are dominated by a single exponential. It is, however, very remark- 
able, that up to t ~ 20 they have rather different slopes, and only beyond that value they all settle 
for the smallest slope. In particular the interpolators with a relative minus sign (derivative interpo- 
lators) show a steeper slope compared to the interpolators with relative plus signs, indicating that 
the former strongly couple to excitations. We stress, however, that also the derivative interpolators 
couple to the ground state, as can be seen in the effective energy plots (r.h.s.): They clearly show 
a second plateau at the ground state energy for large values of t. Interpolator even seems to 
couple to a higher excited state, but settles at the ground state energy beyond t rj 13. The fact, that 
some of our interpolators couple very strongly to excitations shows that they have a large overlap 
with the true physical states. It must, however, be understood, that the situation in 2-d is particularly 
simple due to the lack of angular momentum. All a spatial wave function can do is to have nodes 
in the single spatial direction. Such wave functions are easily obtained by linear combinations of 
our interpolators. 
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Figure 2: Effective energies from the eigenvalues of a 5 x 5 correlation matrix. We compare results from a 
14 x 48 lattice (l.h.s.) to the data from 20 x 48 (r.h.s.), both for g = 0.05, m = 0.05 and N f = 2. 



Let us now come to the analysis of the results from the variational analysis. In Fig. g we 
show the effective energies for the eigenvalues from the generalized eigenvalue problem for a 5 x 5 
correlation matrix built from the interpolators 6)3,05, 08, Gig and O19 for g = 0.05, m = 0.05, 
Nf = 2. We compare two different volumes, 14 x 48 in the l.h.s. plot and 20 x 48 on the r.h.s. 

The plots show long reliable plateaus, in particular for the three lowest states. The two highest 
states on the larger lattice (r.h.s.) show noticeable deviations from a horizontal line, which might 



be an effect of mixing with other states through the correction term in (|1.3|). 

It is interesting to note that except for the lowest lying state (circles) all states show a pro- 
nounced volume dependence and are shifted towards smaller energies for the larger volume. This 
indicates that the higher states might be scattering states, which show, as discussed, a pronounced 
volume dependence. In a future contribution a more detailed analysis of the volume effects will be 
presented, which establishes that indeed the excited states of Fig. ^ are scattering states. 

Here we would like to focus on a different aspect: For the finite volume analysis it is im- 
portant to compare the energies from several different lattice volumes. However, when changing 
the volume, the spectrum is shifted and it is a non-trivial task to identify the individual states, in 
particular since two states might have changed their relative position. For the purpose of matching 



states on different volumes, the eigenvectors of the generalized eigenvalue problem (1.2) are an 
important tool. Since they contain the information on how a state is composed from the individual 
interpolators, they serve as a "fingerprint" for an individual state. 

In Fig. H| we display the entries of the eigenvectors (normalized to 1) for the ground state (l.h.s.) 
and our highest excited state (r.h.s.) as a function of t (16 x 48, g = 0.01, m = 0.05, Nf = 2). 
Also the eigenvector entries form pronounced plateaus and the position of the plateau gives the 
coefficient for the corresponding interpolator in the linear combination which builds up the states. 
These coefficients are surprisingly stable for different volumes and thus are suitable quantities to 
be used in the identification of the states. Furthermore, the coefficients provide physically relevant 
information on how the states are composed: The ground state (l.h.s.) is dominated by the operators 
O3 and O5, whereas the highest excited state (r.h.s.) receives important contributions also from the 
derivative interpolators. Such an analysis allows for a qualitative understanding of the nature of the 
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Figure 3: Entries of the eigenvectors for the ground state (l.h.s) and the fourth excitation (r.h.s) as a function 
of t (16 x 48, g = 0.01, m = 0.05, Nf = 2). The numbers in the legend are labels of the interpolators. 

excited states. In an upcoming paper we will use the methods presented here for a detailed analysis 
of the spectrum of excitations in the lattice Gross Neveu model. 
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